Perturbations of Schwarzschild black holes in laboratories 
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It is well-known that the perturbations of Schwarzschild black holes are governed by a wave 
equation with some effective potential. We consider perturbations of a gas in a tube called the de 
Laval nozzle, which is narrow in the middle and has a sonic point in the throat. By equating the wave 
equation in a de Laval nozzle of an arbitrary form with the wave equation of spin-s perturbations 
of Schwarzschild black holes, we find the exact expression for the form of the de Laval nozzle, for 
which acoustic perturbations of the gas flow corresponds to the general form of perturbations of 
Schwarzschild black holes. This allows observation, in a laboratory, of the acoustic waves, which 
are analogue of damping quasinormal oscillations of Schwarzschild black holes. The found exact 
acoustic analog allows to observe also some other phenomena governed by the wave equation, such 
as the wave scattering and tunneling. 

PACS numbers: 04.30.Nk,04.50.+h 



I. INTRODUCTION 

One of the most essential characteristics of black holes 
are the resonant frequencies of the response to external 
perturbations, called quasinormal frequencies Being 
an analog of normal modes for open systems, they do not 
depend on the way in which the system is perturbed, but 
only on the parameters of a system. Quasinormal modes 
are expected to be observed in the near future with the 
help of a new generation of gravitational antennas. It 
is well-known that non-rotating astrophysical black hole 
can be described by the Schwarzschild solution, implying 
the importance of the quasinormal modes of this back- 
ground 0. 

In addition to the still elusive possibility to observe 
quasinormal modes (QNMs) of black holes with the help 
of a new generation of gravitational antennas, there is 
a window for observation of the acoustic analogue of a 
black hole in laboratories. This is the well-known Unruh 
analogue of black holes [H , which are the apparent hori- 
zons appearing in a fluid with a space-dependent veloc- 
ity, in the presence of sonic points. The supersonic waves 
cannot propagate back beyond the sonic point, mimick- 
ing thereby, the effect of the horizon at sonic points in a 
membrane paradigm. The Unruh discovery stimulated 
active investigation of quasinormal modes of different 
analogue black holes 0]. First of all, the quasinormal 
spectrum of the analogue black holes given by the met- 
rics: 



ds z 



-{l-C/r 2 )dt 2 + (1-C /r 2 )- 1 ^ 2 -2Bdc/)dt+r 2 d(t) 2 



and 



ds 2 = -(1 - r /r 4 )dt 2 + (1 - ro/r 4 )- 1 ^ 2 + rW 
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were considered [6j. These are the two models for rotat- 
ing analogue black holes in "draining bathtub" and for 
a canonical analogue non-rotating black holes. Recently, 
interesting acoustic analogues of brane- world black holes 
were suggested @, [1] ■ 

We can see from the above formulas that these met- 
rics, even being very useful as analogues with apparent 
horizons, do not represent a true solutions of the Ein- 
stein or other gravity dynamical equations. If one had 
a complete analogy with some known solution of Ein- 
stein equations, say, the Schwarzschild solution, he could 
see in the acoustic experiments not only qualitative, but 
also, up to an experimental accuracy, exact numerical co- 
incidence with a prototype characteristics. Namely, for 
quasinormal modes, which are governed by the form of 
the wave equation, this numerical correspondence would 
mean that the effective potential of the perturbations of 
some hydrodynamic system coincides with an effective 
potential of a black hole. Fortunately, recent considera- 
tion of the perturbations of a gas in de Laval nozzle [H 
gives us such an opportunity of finding a system that is 
realted to the same effective potential as a Schwarzschild 
black hole. 

The canonical de Laval nozzle is a convergent-divergent 
tube, narrow in the middle. It allows to accelerate the 
gas until the sonic speed in its throat, reaching super- 
sonic speeds after passing the throat. The perturbations 
of the gas in de Laval nozzle can be considered as one- 
dimensional if the section of the nozzle does not change 
too quickly. Here we show that the corresponding effec- 
tive potential of perturbations in a canonical de Laval 
nozzle can be made to be equal to the potential for per- 
turbations of Schwarzschild black holes, if choosing some 
specific form of the nozzle. In addition, we suggest an- 
other, approximate, way to get quasinormal modes of 
Schwarzschild black holes in de Laval nozzle. For this, 
one needs to mimic the form of the effective potential for 
Schwarzschild metric with the help of a de Laval nozzle 
of a simple form suggested in [1] . 

The paper is organized as follows: Sec HT1 gives all ba- 



2 



sic equations of the one-dimensional motion in de Laval 
nozzle we shall use. Sec [Hi] is devoted to reproducing 
the exact expression for the form of the nozzle which 
corresponds to the potential of the Schwarzschild black 
holes. In the discussions, we sketch the open questions 
and possible generalizations of the suggested technique. 



where M{x) is the Mach number which, by the def- 
inition, is the current flow speed divided by the sound 
speed. In our notations M = v/c s . The function H u 
represents small perturbations of gas flow, 



H u (x)=9 1/2 



dt e l ^~ f ^cj){t,x), (12) 



II. CALCULATION OF THE NOZZLE CROSS 
SECTION A IN TERMS OF g 

We assume that a gas in the nozzle can be described 
by equations of motion for perfect fluid and that the flow 
is quasi-one-dimensional: 



d t (pA) + d x (pvA) = , 
dt(pvA) + d x [( P v 2 + P )A] =0, 
d t (eA)+d x [{e + p)vA] =0, 



(1) 
(2) 
(3) 



Here p is the density, v is the fluid velocity, p is the 
pressure, A is the cross section of the nozzle, and 
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P 



1 



(4) 



is the energy density. The heat capacity ratio is 7 = 
1 + 2/n = 7/5 = 1.4 for di- atomic molecules of air 
(n = 5). We shall assume that the flow has no entropy 
discontinuity, then the fluid is isentropic 

p oc p 1 . (5) 

Instead of Eq. ([2]), we can use Euler's equation 

p(d t + vd x )v = -d x p , (6) 

is reduced to the Bernoulli's 



For isentropic fluid Eq 
equation 

1 



d t $ + -(d x <py + h( P ) = o , 



(7) 



where h(p) = J p~ 1 dp is the specific enthalpy and $ = 
J v dx is the velocity potential. 

According to [J] , the perturbation equations in such a 
nozzle can be reduced to: 



dx* 



V(x*) 

CsO ' 



0, 



V(x*) = - 



a d 2 g 



2 dx* 



4 \dx*J 



(8) 
(9) 
(10) 



Here c s q is the stagnation sound speed, and x* is an 
acoustic analogue of the tortoise coordinate which sat- 
isfies x*(x = +00) = +00, x*(x = 0) = —00, namely, 



CsO 



dx 



c s {\-M( x yy 



(11) 



9 = —, 



(13) 
(14) 



Here, according to [J], the small perturbations are defined 
as follows 

p = p + Sp, p^\Sp\, (15) 
$ = $ + 0, (16) 

Our starting point is the calculation of the configura- 
tion of de Laval nozzle, i. e. its cross section as a func- 
tion of the transversal nozzle coordinate. Since we know 
the effective potential, we can calculate in some way the 
function g(x). By definition Eq. (15) of [3| 

<r pA 

Cs \fjpfp 

Taking (5) into account we find 

P A 



g oc 



p(7-l)/2' 



(17) 



We can choose dimensionless quantities for p(x) and 
A[x) by measuring them in units of po and A* respec- 
tively @. Then equation (5.3) of [5[ reads 



A 



^cx (l-p^) 1 ' 2 p 



(18) 



Since ([5]) is invariant with respect to re-scaling of g, we 
can fix the coefficients in (fl~7| and (fT8|) arbitrarily: 



P A 



2 j0 ^- 1 )/ 2 ' 



A 



1 = (l - P^) 1/2 P. (19) 



We find 



,(l-7)/2 

T72 



P 



(1-7) 



.1/2 



(20) 



2(1-/^-1))^ 2(p( 1 -^)-l) 
Hence it follows that 

P ^ = 2. 9 2 (1 ± yr^F 1 ) • (21) 

The sign should be chosen in order that p be a 
monotonous function with respect to the transverse co- 
ordinate. As we will show later, the function g for the 
Schwarzschild black hole can be chosen also monotonous 
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in the R region, finite at the horizon and infinite at the 
spatial infinity. Therefore, we choose the minus sign, 



1 -' = 2 fl a l-vW 



(22) 



Note that g must always be larger than unity in our con- 
sideration. 

In our notations, the Mach number is connected with 
p as 



We find 



P 



M 



(23) 



1 (P 1 - 1 ~ 1) 



7- 



T ( 2 . 9 2 (l-v^F^)-l) 



(24) 



Since M = 1 at the event horizon, g must be finite there, 
and 



7+1 



= -=>!. 



(25) 



This requirement fixes both constants of integration. 

Substituting (j2"3")l in (|19p we find the cross-section area 
as a function of g: 



A 



V2 [2g 2 (l - 



1/(7-1) 



(26) 



III. THE DE LAVAL NOZZLE FOR THE 
SCHWARZSCHILD BLACK HOLE 

Our starting point is the calculation of the configura- 
tion of the de Laval nozzle, i. e. its cross section as a 
function of the transversal nozzle coordinate. Since we 
know the effective potential, we can calculate in some 
way the function g(x). 

After separation of the angular and time variables, 
scalar field perturbations in the Schwarzschild back- 
ground, can be reduced to the wave-like equation 



(27) 



— +u, z -V(r*))V(r*) = 0, 



where putting the event horizon to be unity we find, 

1 . , dr 



f(r) 



1 



V(r) = f(r) 



dr* 



1) 



f(r) 



1-s 2 



(28) 



To find g that produces the same potential we identify 
the "tortoise" coordinates of the black hole solution and 
of the laval nozzle 



dr* = dx* 



c s0 dx p 0--yV 2 dx 



1 - M 2 



2g 2 (l-^/T^g^)dx 



1 



(29) 



Here x is the real coordinate along de Laval nozzle. Then 
we can find the equation for g(r), 



f(r)f(r)g'(r)+f(r) 2 g"(r) f(r) 2 g'(r) 



Mr) 



4.9 (r) 5 



V(r). (30) 



This implies that the from of de Laval nozzle is parame- 
terized by the parameter r and x*(r) — r*(r). Note that 
as we chose the radius of the event horizon to be unity, 
the nozzle coordinate x is measured in the units of the 
radius of the event horizon. 

The general solution of the equation (|30j) contains two 
arbitrary constants. They can be fixed in a unique way 
by the condition (|2"5|) . Namely, the requirement that the 
solution must be finite at r = 1 fixes one of the constant. 
Then the other constant re-scales the solution of (|30|) . 
and must be fixed by its value at r = 1. Finally, the 
solution of (|3"0|) for arbitrary I and s, that satisfies ([23)1 
is given by the following formula: 

9(r) = 



7 + 1 



(-l) n+s (£ + n)\ 



„n+l 



272^7^1 ^ V (« + s)l(n - s)l(t - n)\ 

= 7 + 1 r 2*+2 x 

r(l + l + s) 2 F 1 (s-£,s + £+l,l + 2s,r) 



(31) 



r(i + 1 - s)r(i + 2s) 

One can easily check that the above solution indeed sat- 
isfies the equation (|25|) . for any fixed £ and s. 

From (f2"9")l we find the dependance of the transversal 
nozzle coordinate x on the parameter r: 



(7 + 1 - 4g(r) 2 (l - Vl-ff(r)-*)) dr 
I /(r)( 7 - l)^2g{rf (l - ^T^gW^) 



(32) 



The integration constant is chosen in order to x be zero 
at the sonic point. 

Now we are in position to find the required form of 
de Laval nozzle, i.e. to find its cross-section A{x). We 
just need to replace g(r) given in (f3"Tj) in (|2"6"|) and go 
over to the transverse nozzle coordinate x. The function 
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Figure 1: The form of de Laval nozzle and the effective potential for s = £ — 0. 




Figure 2: The form of de Laval nozzle and the effective potential for s = £ = 1. 
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A(x) is shown in Figs. 1-4. Note that the canoni- 
cal de Laval nozzle is diverging at the end of the flow 
trajectory, so that A x=OD — oo (see pages 53 and 124 
in [a]). Indeed, our formula (j3Tj) implies divergence at 
least as ~ r 2 . The diverging of the nozzle nevertheless 
does not give any going beyond the one-dimensional rep- 
resentation of the motion, because the function y A(x) 
is measured in units of black hole mass, i.e. one can 
" pull" the nozzle along the transverse coordinate x in or- 
der to make the area of the nozzle change as slowly as one 
wishes. Such a "pulling" simply means that we are get- 
ting the correspondence with a black hole of larger mass. 
Since the quasinormal modes are inversely proportional 
to the mass of the black hole, this means just some deter- 
mined multiplication by a coefficient when coming from 
the frequencies observed in experiment to the QNMs of 
a black hole. 

As can be seen from |9[ , the values of the quasinormal 
modes arc determined by the behavior of the effective po- 
tential in some region near black hole. The form of the 
effective potential (and thereby of de Laval nozzle) far 
from black hole is less significant for the QNMs problem. 
Therefore we expect that such experimental phenomena 
as surface friction and reflection of waves from bound- 
aries will not have considerable influence on the observed 
picture. 

Now we discuss the isospectrality and the effective po- 
tential for the polar and axial gravitational perturba- 
tions. We consider the effective potential for the grav- 
itational perturbation of the polar type, 



V(r) = f(r) 



9(1 + Xr) + A 3 r 3 + A 2 r 2 (3 + 2r) 
r 3 (3 + A) 2 • 



where A = (£ + 2)(i-l). 

For the above polar type gravitational perturbations 
we can also obtain the exact solution for the function 
g(r), i.e. the form of de Laval nozzle. The function g(r) 
is given in the following table for £ — 2, 3, 4, and 5 and 
in the formula 



9(rJ) 



7 + 1 



r 2 p(r, £) 2 



2^Vt^T (3 + (^ + 2)(^-l)r) 2 



(34) 



i 


p(r,£) 


2 


3 - 6r 2 - 4r 3 




3 


3 - 30r 2 - 20r 3 + 


60r 4 


4 


3 - 90r 2 - 60r 3 + 


630r 4 - 504r 5 


5 


3 - 210r 2 - 140r 3 


+ 3570r 4 - 6552r 5 + 3360r 6 



(33) 



Apparently there exists a solution for general t, yet 
probably quite cumbersome. Analysis of the Figs. 3-4 
shows that the form of the nozzles for modeling polar and 
axial gravitational perturbations are almost the same. 
The difference cannot be seen explicitly, although not 
vanishing. As the effective potentials for axial and polar 
types also differ only slightly, this may mean that the 
forms of the nozzles also differ only slightly. 



IV. DISCUSSION 

The suggested solution of the inverse problem for the 
correspondence of the form of de Laval nozzle to the 
general form of perturbations of the Schwarzschild black 
holes (i.e. for perturbations with spin s and multipole £) 
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Figure 3: The form of de Laval nozzle and the effective potential for s = £ — 2. 
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Figure 4: The form of de Laval nozzle and the effective potential for the polar gravitational perturbation 1 = 2. 
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can be generalized in many ways. First of all, it would be 
very interesting to consider the massive vector [TlT | and 
scalar (lpj field perturbations, because of quite unusual 
behavior of massive perturbations. Thus the so-called 
quasi-resonances, infinitely long lived modes, for massive 
scalar field [11] could be observed in a de Laval nozzle of 
some form as almost non-damping sound waves. These 
waves would certainly be reflected from the boundary of 
the nozzle and, thereby, would break the quasinormal 
mode boundary conditions. Even though in a real exper- 
iment one cannot obtain perfect QNM boundary condi- 
tion (QNM b.c), a considerable deviation from QNM b.c. 
should be observed when modeling the quasi-resonances. 

Another possible generalization is to consider 
more general black hole backgrounds: Reissner- 
Nordstrom, Schwarzschild-de Sitter, or higher dimen- 
sional Schwarzschild black holes [13[ with charge, A-term 
and Gauss-Bonnet-term [l5| . including the brane- world 
black holes [TiJ. The Gauss-Bonnet black hole ringing 
[l5| would be especially interesting to model in a nozzle 
because there exists the instability in some region of 
values of the black hole parameters [15| . In the approach 
considered in this paper we are limited only by a spher- 
ical symmetry, i.e. by lu independence of the effective 
potential. In addition, one could consider the flow of gas 
with a time dependent initial speed at the compressor, 
which probably could model the perturbations of the 



Vaidya evaporating black holes [l6j], 01- We cannot be 
sure, that in all these cases the differential equation for 
the form of de Laval nozzle will be exactly integrable, 
yet one can always find some numerical solution. We 
believe that further research will solve these interesting 
problems. 

It should be recalled also that the precise acoustic anal- 
ogy is only established for a scalar field. To be able to 
reproduce the potential V(r) for fields of different spins 
certainly does not mean that one can reproduce all the 
characteristics of those equations in an acoustic model. 

Finally, let us note that the obtained acoustic ana- 
logue for the perturbations of the Schwarzschild black 
holes is not limited by quasinormal mode problems only, 
but allow general investigation of propagation of classical 
and quantum fields, including such processes as scatter- 
ing and tunnelling of waves and particles. 
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